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ABSTRACT 

The Centaurs are recent escapees from the Kuiper belt that are destined 
either to meet fiery oblivion in the hot inner regions of the Solar system or to 
be ejected to the interstellar medium by gravitational scattering from the giant 
planets. Dynamically evolved Centaurs, when captured by Jupiter and close 
enough to the Sun for near-surface water ice to sublimate, are conventionally 
labeled as "short-period" (specifically, Jupiter-family) comets. Remarkably, some 
Centaurs show comet-like activity even when far beyond the orbit of Jupiter, 
suggesting mass-loss driven by a process other than the sublimation of water ice. 
We observed a sample of 23 Centaurs and found nine to be active, with mass-loss 
rates measured from several kg s~^ to several tonnes s~^. Considered as a group, 
we find that the "active Centaurs" in our sample have perihelia smaller than 
the inactive Centaurs (median 5.9 AU vs. 8.7 AU), and smaller than the median 
perihelion distance computed for all known Centaurs (12.4 AU). This suggests 
that their activity is thermally driven. We consider several possibilities for the 
origin of the mass-loss from the active Centaurs. Most arc too cold for activity 
at the observed levels to originate via the sublimation of crystalline water ice. 
Solid carbon monoxide and carbon dioxide have the opposite problem: they are so 
volatile that they should drive activity in Centaurs at much larger distances than 
observed. We consider the possibihty that activity in the Centaurs is triggered 
by the conversion of amorphous ice into the crystalline form accompanied by the 
release of trapped gases, including carbon monoxide. By imposing the condition 
that crystallization should occur when the crystallization time is shorter than the 
orbital period we find a qualitative match to the perihelion distribution of the 
active Centaurs and conclude that the data are consistent with the hypothesis 
that the Centaurs contain amorphous ice. 

Subject headings: comets, Kuiper Belt, cometsrgeneral, asteroids: general 
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1. Introduction 

The Centaurs experience repeated close encounters with the giant planets and, as a 
result, follow orbits which are both chaotic and dynamically short-lived compared to the 
age of the Solar system. Numerical simulations give half-lives from 0.5 Myr to 30 Myr, 
with average values increasing with semimajor axis and perihelion distance (Horner et al. 
2004) and a population median hfetime near 10 Myr (Tiscareno and Malhotra 2003). The 
short dynamical lifetimes suggest that the Centaurs must be resupplied from a more stable 
source region elsewhere, if their population is to remain in steady-state. This source region is 
most likely the Kuiper belt, from which objects drift first into Neptune-crossing trajectories 
and then are scattered between the giant planets (see Morbidelli 2008 for a recent review 
of Kuiper belt dynamics). The Kuiper belt contains several dynamically distinct classes 
of body, and it is not clear precisely from where in the belt the Centaurs are derived. The 
source is popularly identified with the so-called scattered Kuiper belt objects but other trans- 
Neptunian populations, or even completely different source regions, may be implicated (c.f. 
Volk and Malhotra 2008). The Centaurs suffer several fates. A few percent of the Centaurs 
impact the giant planets, or are broken apart by gravitational stresses exerted upon them by 
near-miss interactions with the giants (the famous comet D /Shoemaker-Levy 9 is the best 
example of a body suffering both of these fates). About two thirds are ultimately ejected 
from the Solar system into interstellar space (Tiscareno and Malhotra 2003). The remaining 
one third become temporarily trapped in small, near-Sun orbits, where they sublimate and 
are conventionally labeled as "Jupiter family comets" . The Centaurs, then, are scientifically 
interesting as a kind of intermediate population, halfway between the distant, frigid and 
relatively stable bodies in the Kuiper belt and the rapidly sublimating, disintegrating comets 
of the hot inner regions (see Jewitt 2004 for a broad account of the interrelations between 
these objects). 

Surprisingly, there is no uniformly accepted definition of the Centaurs. Some dynami- 
cists classify all objects with perihelia between the orbits of Jupiter and Neptune as Centaurs, 
regardless of their semimajor axes. This definition is reasonable, given that the dynamics 
are dominated by perihelic interactions with the planets but some qualifications are needed, 
since many resonant Kuiper belt objects (like 134340 Pluto), have perihelia inside Neptune's 
orbit but are not usefully considered as Centaurs. A detailed classification system has been 
proposed by Horner et al. (2004) in which the Centaurs are labeled by the planets with which 
they most strongly interact. We here adopt a simpler operational definition, modified from 
Jewitt and Kalas (1998), and suggested by the prototype Centaur, (2060) Chiron (Kowal 
et al. 1979), whose dynamics are dominated by perihelic and aphelic interactions with giant 
planets. We define the Centaurs as bodies whose orbits meet the following conditions: 
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• The perihelion distance, q, and the semimajor axis, a, satisfy aj < q < and 
aj < a < Oat, respectively, where aj = 5.2 AU is the semimajor axis of Jupiter and ojv 
= 30.0 AU is the semimajor axis of Neptune. 

• The orbit is not in 1:1 mean-motion resonance with any planet. 

This definition has the practical advantage of simplicity; it efficiently isolates a sub-set 
of the Solar system small-body population whose members are short-lived by virtue of their 
gravitational interactions with the giant planets. At the time of writing there are, by the 
above definition, 92 known Centaurs. 

In this paper, we focus attention on Centaurs which display comet-like mass-loss. Our 
purpose is to identify common properties of the activity and to seek evidence concerning 
its origin. The active Centaurs have, so far, escaped systematic study as a group. One 
likely reason is that the cometary classification afforded to such objects makes no distinction 
between them and other, long-studied comets in the Jupiter family, Halley family and Long- 
Period populations. We seek answers to such questions as 1) How do the properties of active 
Centaurs compare with those of other comets and with Centaurs as a whole? 2) Is mass-loss 
implicated in the reportedly bimodal distribution of the optical colors of the Centaurs? 3) 
What drives the mass-loss from the active Centaurs? 



2. Observations 

The imaging observations reported here were taken at the University of Hawaii 2.2-m 
telescope (UH 2.2-m) and the Keck 10-m telescope, both located atop Mauna Kea, Hawan. 
At the UH 2.2-m we employed a Tektronix 2048x2048 pixel charge-coupled device (CCD) 
camera having an image scale of 0.219 arcsec (pixel)~^ when mounted at the f/10 Cassegrain 
focus. At the Keck, we employed the Low Resolution Imaging Spectrograph (LRIS) which, in 
its modern form, employs independent red and blue channels separated by a dichroic beam- 
sphtter (Oke et al. 1995). A 5000 A dichroic was used. The red-side detector is a front-side 
illuminated Tektronix 2048x2048 pixel CCD camera, with image scale 0.210 arcsec (pixel)"^ 
and two read-out amplifiers. The blue-side used two Marconi (E2V) CCDs, each 2048x4096 
pixels in size with two readout amplifiers per chip, and with a scale of 0.135 arcsec (pixel)"^. 
On the blue side we used a B filter with effective wavelength Ag = 4370A and full-width at 
half maximum (FWHM) of AA = 900A. On the red side, the filters employed were V (Ag = 
5473A, AA = 950A), R (Ae = 68OOA, AA = 1270A) and I (Ae = 8330A, AA = 3130A). 

At both telescopes, we constructed fiats by taking dithered images of the twilight sky 
and then median combining these images, after subtracting a bias frame and scaling them 
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to a common value, to eliminate stars. For comparison, we also constructed "data fiats" 
from the median combination of images of astronomical targets taken through the night, 
relying on the lack of spatial correlation between images to eliminate field stars. We found 
that the data-flats were superior to the twilight-sky flats in the B and I filters and we used 
them exclusively. In V and R the two types of fiat were comparable in quality. Photometric 
calibration of the flattened images was obtained from measurements of Landolt standard 
stars (Landolt 1992). We used standard stars with optical colors similar to, or slightly 
redder than, the Sun, in order to minimize color terms. The faintest Landolt stars were 
used in order to minimize any possible shutter errors, and we experimented with different 
exposures to confirm that the shutter timing response was linear at the 1% (0.01 mag.) 
level. Lastly, we avoided standard stars with listed photometric uncertainties greater than 
0.01 mag. 

We observed a total of 23 Centaurs. A journal of observations is given in Table [T] and 
the main orbital parameters of these Centaurs are listed in Table [2j Not all objects listed in 
Table [T] appear in subsequent Tables because of field-star contamination and other problems 
that caused us to reject some data from further consideration. The estimated semimajor axes 
of two of the objects in our sample (namely, (148975) 2001 XA255 with a = 30.183 AU and 
2002 PN34 with a = 30.781) have been revised following the acquisition of observations, such 
that they (slightly) exceed the semimajor axis of Neptune (30.070 AU). We have elected to 
keep them in the sample, while noting that none of the conclusions drawn would be different 
if we instead chose to exclude them. 

The next section gives brief introductions to the active Centaurs considered in this work. 

2.1. 29P/Schwassmann-Wachmann 1 

This object is labeled, and traditionally regarded, as a comet but it also meets the 
dynamical definition of a Centaur and so it is appropriate to consider it here. Our optical 
observations (Table [3]) show 29P in an active state, with a visible extent much larger than 
100 arcsec (Figure [T]). The colors are slightly redder than the Sun and essentially constant 
with respect to aperture size, except for a marginally significant trend towards bluer colors 
at larger radii. Earlier observers have reported that the broadband colors of 29P measured 
within a very large (88 arcsec diameter) aperture vary with the level of activity, becoming 
redder as the comet brightens (Kiselev and Chernova 1979). 

Since discovery in 1927, 29P has shown repeated photometric outbursts superimposed 
upon a platform of constant activity (e.g. Roemer 1958, Jewitt 1990, Trigo-Rodriguez et al. 
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2008). Spectroscopic measurements show that the activity is driven by the subhmation of 
carbon monoxide (CO) ice (Senay and Jewitt 1994; Crovisier et ah 1995), with CO mass- 
loss rates of 1000 to 2000 kg s'^ Other species have been reported when the comet is at its 
brightest, including CO^ and CN (Cochran and Cochran 1991). 

Thermal measurements give an upper limit to the nucleus radius of 27±5 km and a 
lower limit to the geometric albedo of pr = 0.025±0.01 (Stansberry et al. 2004). (These 
values are limits because of the substantial contamination of the infrared data by coma, and 
the need to make a model-dependent correction for coma.) 

2.2. 39P/Oterma 

The orbital parameters of this Centaur evolve substantially on timescales of decades 
owing to close interactions with Jupiter. For example, the perihelion was raised from 3.39 
AU to 5.47 AU by a close approach (0.1 AU) to Jupiter in 1963, after which 39P was not 
re-observed until 1998. In our data from the UH 2.2-m telescope taken 2002 Sep 07, 39P was 
found to show a weak coma, barely visible in Figure [2j The object is slightly redder than 
the Sun, with no compelling evidence for a color gradient within the apertures employed for 
photometry (Table [s]). 

2.3. 165P/Linear 

Coma was reported around this object (formerly known as 2000 B4) in observations 
taken with a 0.65-meter diameter telescope on UT 2000 Feb 10 (Kusnirak and Balam 2000) 
but no coma has been reported since this initial report. Our own attempts to observe the 
object were made on UT 2002 Jan 09 using the UH 2.2-m telescope and on UT 2003 Jan 08, 
using the Keck 10-m telescope. The object was not detected down to approximate limiting 
magnitudes mn = 23 in the former case and 24.5 in the latter. We conclude that either the 
object was very faint or the ephemeris position was in error by more than the field radius 
(about 3.5 arcminutes). Object 165P/Linear will not be further discussed in this paper. 

2.4. 166P/2001 T4 

Our imaging data show a clear coma at all epochs (Figure [3]) and the BVRI color data 
show that 166P is optically very red. For example, in 2001 Nov (Table |5]) we find B-V 
= 0.87, V-R = 0.69, R-I = 0.78 within a 5 pixel (1.1 arcsec) radius aperture. Bauer et 
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al. (2003) reported 166P to be one of the reddest objects in the Solar system, with V-R = 
0.95±0.04 and R-I = 0.69±0.07 measured in 2001 Oct and 2002 Jan within a 1.5 arcsec radius 
aperture. The reason for the difference is not clear: we used the same telescope and filters and 
our measurements were obtained within a month of Bauer et al.'s first observations. These 
authors also reported a radial color gradient in 166P, with the image core being substantially 
redder than the coma. Our data confirm the existence of such a radial color gradient (Table 

2.5. 167P/2004 PY42 

This Centaur was observed on UT 2004 Oct 10 at the Keck 10-m telescope under 
photometric conditions. A weak coma was apparent in the raw data, and can be seen in 
Figure |4] (note that the diffuse object to the north east of 167P in the figure is a background 
galaxy). The coma is faint compared to the nucleus. 

2.6. C/2001 MIO (Neat) 

This Centaur was observed on UT 2002 Sep 08 at the UH 2.2-m telescope and on UT 
2004 Oct 10 using the Keck 10-m telescope (Figure |5]). As a result of instrumental problems, 
observations on the former date were obtained only in the B and R filters, while on the latter 
date full BVRI colors were secured. The object faded by ~3.5 magnitudes between 2002 and 
2004, as the heliocentric distance increased from 6.056 AU to 9.145 AU. Only ~2 magnitudes 
of this fading can be attributed to geometric effects, the remaining ~1.5 mag. reflects a real 
decrease in the coma cross-section of C/2001 MIO. However, the B — R color indices at the 
two epochs of observation, 1.22±0.08 and 1.12±0.09, are completely statistically consistent 
within the measurement uncertainties. Since the coma contribution to the signal in the UT 
2004 Oct 10 data was reduced relative to that from the nucleus by a factor ~4, this suggests 
either that the nucleus contributes a negligible fraction of the total cross-section even at 9 
AU or that the nucleus and the coma have the same B — R color. 

2.7. P/2004 Al (Loneos) 

This object passed within 0.03 AU of Saturn in 1992 causing a sudden decrease in 
the perihelion distance from ~10 AU to ~5 AU (Hahn et al. 2006). We observed it post- 
perihelion on UT 2006 July 01 at i? = 6.079 AU, when it displayed a prominent coma (Figure 
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[6]). The absolute magnitude we derived, 1,0) > 10.92, is consistent with the faintest 

values reported by Mazzotta et al. (2006). 

2.8. 174P/Echeclus (60558) 

The object (formerly asteroid 2000 EC98) was discovered and measured as an inert 
Centaur (Rousselot et al. 2005) but promptly erupted into activity in late 2005 (Choi and 
Weissman 2006). We imaged the object using Keck on 2006 February 25 (Figure |8]). The 
coma showed a complicated post-outburst structure, visible up to 65 arcsec from the nucleus. 
We digitally removed background stars from the images and determined magnitudes within 
the standard apertures. The magnitude of the nucleus corrected to unit heliocentric and 
geocentric distances and to zero degrees phase angle, mR(l, 1, 0) > 8.7 (Tablejsj), is consistent 
with the pre-activity value, Hji = 9.03±0.01, measured by Rousselot et al. Only R filter color 
measurements were possible at Keck due to time constraints (Table [s]). However, spatially 
resolved colors were determined on UT 2006 Feb 24 by Bauer et al. (2008). They found a 
coma slightly redder than sunlight with no substantial color gradients in the 5 arcsec to 30 
arcsec projected radius range. 

The nature of this object remains a mystery. The motion of the active secondary compo- 
nent suggests that it is a fragment ejected from the parent body, but neither radiation forces 
nor non-gravitational forces due to asymmetric outgassing seem capable of producing the 
large relative velocity of separation (Weissman et al. 2006). Thermal infrared measurements 
show that the coma particles are large and emitted continuously, not impulsively (Bauer et 
al. 2008). 

2.9. (5145) Pholus 

(5145) Pholus (formerly 1992 AD) is well known for its ultra-red optical reflection spec- 
trum and presumed organic composition. No coma has ever been reported around this object 
and no coma was detected in our Keck images taken UT 2004 Feb 17 (Table [s]). 

2.10. (63252) 2001 BL41 

This object was observed as a point source on UT 2002 Feb 22 and 23 at the UH 2.2-m 
telescope with very similar apparent magnitudes on the two nights (Table [3]). 
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2.11. 2003 QD112 

We observed 2003 QD112 at the Keck telescope on UT 2004 Oct 10 when at 12.65 AU 
hehocentric distance and discovered a coma about this object (Figure [?]). Only i?-band data 
were secured (Table |3|). We reimaged 2003 QD112 at the Keck on UT 2007 Oct 12. A coma 
was again evident in a blinked pair of 300-sec R-band integrations, but the comet appeared 
in projection near to a much brighter field galaxy, making accurate photometry impossible 
on this date. 

2.12. P/2005 S2 (Skiff) 

We attempted to image P/2005 S2 on UT 2007 Sep 19 using the UH 2.2-m telescope. 
The object was not detected in 300-s i?-band exposures with a limiting magnitude near R 
= 23.5, leading us to conclude either that the brightness was less than expected or that the 
ephemeris was substantially in error. Subsequently, we detected P/2005 S2 as a point-source 
object near mag nitude 23.0 (Table |3|) using the Keck 10-m on UT 2007 Oct 12. 

2.13. P/2005 T3 (Read) 

This Centaur was not detected on either of two attempts (see Table [T]) with two different 
telescopes. On UT 2007 Sep 17 the object was projected against a background H2 nebula 
which may have impeded its detection. On UT 2007 Oct 12 a bright background star filled 
the field of view with scattered light. 

3. Observational Results 

3.1. Coma Search 

We first visually searched each Centaur for coma by eye, comparing images taken at dif- 
ferent times in order to distinguish faint near-nucleus structure from occasional background 
objects. We also used multi-aperture photometry to measure, or set limits to, light from 
near-nucleus coma. Given a point-source nucleus embedded in an extended coma, the in- 
ner aperture of radius 0i would sample the Centaur nucleus combined with an uncertain 
contribution from near-nucleus coma while the annulus between 0i and an outer aperture 
of radius 02 would provide a measure of pure coma. Then, the apparent magnitude of the 
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coma in the annulus between 0i and (f)2 is given by 

rrid = -2.5 log^o (10-°-^'"^ - lO""-^™^) (1) 

in which mi and 1712 are the magnitudes within apertures of radii (pi and (j)2- In practice, 
a small correction is needed because the point-spread function (PSF) of the image is non- 
zero in the annulus between 0i and 02. We measured the magnitude of this correction to 
the brightness in the coma annulus using the PSFs of field stars in the same images as the 
Centaurs. 

Based on trial and error, we selected aperture radii 0i = 2.2 arcsec and 02 = 3.3 arcsec 
as standard for our coma-search measurements. Smaller values of 0i proved overly sensitive 
to the steep surface brightness gradients found near the image core and to small differences 
in the PSF caused by relative motion between the Centaurs and the field stars used for 
reference. Conversely, substantially larger apertures, 02? were more frequently contaminated 
by faint background sources and more subject to photometric errors from the uncertainty in 
the brightness of the night sky. With this choice of 0i and 02, we found that the wings of 
the PSF between 0i < < 02 typically contained only ~1% to 2% of the light with < 0i. 

The resulting coma annulus magnitudes, m^, are listed in column 5 of Table [3j Where 
no significant coma is detected, we have listed a 3a limit to computed from the scatter 
in multi-annulus measurements of ~5 to ~10 field stars. Coma is detected in nine of the 
23 Centaurs observed. In eight of these nine cases, coma that was detected through multi- 
aperture photometry and Equation ([T| was also visible to the eye in highly stretched images 
blinked rapidly on a computer. Only in Centaur 2006 SX368 did we find evidence for coma 
from the photometry that could not be corroborated visually. 

3.2. Absolute Magnitudes and Sizes 

The magnitudes of the Centaurs are corrected to Rau = = 1 and phase angle a = 
using 

mR{l, 1, 0) = - 5 logioiRau^au) - $(«) (2) 

where Rau and Aau are the heliocentric and geocentric distances in AU and is the 

so-called phase function correction at phase angle a. The phase angles, a, at which the 
Centaurs were observed span the range from < a < 10° (Table fl]). To first order, the 
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phase function over this modest range of angles can be represented as having a hnear (in 
magnitude space) dependence on a, given by $(a) = l0(-0-4/3"), where (3 [ mag. deg 1 is 
a constant. However, at small phase angles, the phase functions of many objects display a 
narrow-angle peak (the "opposition-surge") while at larger phase angles the phase function 
often appears measurably non-linear. The phase functions of the active Centaurs are further 
uncertain, because the angle dependence of the scattering from the dust particles in the coma 
is not necessarily like that of the solid surface of the nucleus from which those particles were 
ejected. Here, we used the "HG" formalism of Bowell et al. (1989) with scattering parameter 
G = 0.15. This choice gives phase darkening similar to the constant (3 = 0.1 mag. deg~^ 
case (e.g. the Bowell formalism gives phase darkening at a = 10° to 35% of the opposition 
value while /3 = 0.1 mag. deg~^ gives a 40% darkening). 

The absolute magnitudes, mji{l, 1, 0), are listed in Table [sj The Table also lists effective 
object radii, [km], computed under the assumption of geometric albedo pr = 0.1. The 
statistical uncertainty on 771/2(1,1,0) is the same as on m^j but the actual uncertainty is 
larger, owing to uncertainty introduced by the phase angle correction. The uncertainty on 

is potentially large and unknown, pending future determinations of the albedos and phase 
functions of the observed objects. We include 7/1^(1, 1, 0) and Ve here mainly to provide an 
approximate indication of the relative sizes of the observed objects. Note that entries for 
active Centaurs in Table |3] are presented as lower and upper limits on mji{l, 1,0) and r^, 
respectively, with no attempt to remove coma from the central aperture using profile-fitting 
or other techniques. The smallest observed Centaurs, P/2005 S2 (rg = 2 km), 2000 GM137 
(re = 3 km) and 39P (rg < 3 km) have sizes consistent with the nuclei of better-observed 
Jupiter family comets. However, the median radius of the observed (inactive) Centaurs is 

= 15 km (Table [3]), about an order of magnitude larger than the characteristic radii of 
well-measured Jupiter family comet nuclei (Lamy et al. 2004) and an order of magnitude 
smaller than the characteristic radii of Kuiper belt objects for which reliable physical data 
exist. These differences simply reflect the relative mean distances of the different types of 
object and the limitations of magnitude-limited photometry. 



3.3. Colors 

Like the Kuiper belt objects, the Centaurs are remarkable for their wide range of optical 
colors, ranging from nearly neutral to very red. In fact, two of the most famous Centaurs 
span the full range of colors present in the outer Solar system: (2060) Chiron is nearly 
neutral (Meech and Belton 1990, Luu and Jewitt 1990) and shows weak bands of water ice 
(Foster et al. 1999, Luu et al. 2000) while (5145) Pholus is very red (Mueller et al. 1992), 
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perhaps because of the presence of irradiated complex organic compounds (Luu et al. 1994, 
Cruikshank et al. 1998). The fact that Chiron is an active Centaur while Pholus is not has 
lead several investigators to speculate that outgassing activity might be responsible for the 
neutral colors of some Centaurs. For example, outgassing could blanket the surface of a body 
with neutral-colored fine particles excavated from beneath the few meter-thick surface layer 
in which radiation damage is strong (Delsanti et al. 2004). There is some evidence that the 
Centaurs as a whole occupy a bimodal distribution of optical colors (Peixinho et al. 2003). 
If so, the bimodality could tell us which Centaurs have been recently active. However, other 
interpretations are possible. Perhaps activity does not affect the surface colors. The color 
differences could reflect different intrinsic compositions, possibly related to different Centaur 
source regions. 

Figure [9] shows the distribution of the Centaurs in the B-V vs. V-R color plane. The data 
on inactive Centaurs (red circles) are taken from the compilation of Peixinho et al. (2003), 
while data on active Centaurs (blue circles with error bars) are from Table |6} Measurements 
of two Centaurs (Chiron and 60558) appearing in both the Peixinho et al. compilation and 
the present work are joined by lines. For clarity, only the uncertainties on the photometry 
of the active Centaurs are shown. The Figure shows the "red" and "blue" color clumps 
occupied by the Centaurs in the Peixinho sample. Six of the active Centaurs for which we 
possess quality B-V and V-R colors fall within the "blue" clump, while the seventh, 166P, 
is red in V-R but not in B-V, and falls in neither the blue clump nor the red. 

An observation relevant to the more neutral colors of the active Centaurs is provided 



by Figure 10, which shows the B — I color index measured as a function of the projected 
angular distance from the photocenter of the five active Centaurs. Only these five objects 
have adequate spatial data and signal-to-noise ratio to justify an examination of the radial 
color gradient. The wide range in the B — I colors of the active Centaurs is apparent from 
this Figure, from the very red 166P/2001 T4 to nearly neutral C/2001 MIO (the Solar color 
index is {B — I)q = 1.38). The main trend apparent in Figure [l0| is the tendency for the 
B — I color index to decrease with angular radius, meaning that the comae of the active 
Centaurs are typically bluer than their centers. We used the profiles of field stars to confirm 



that the color gradients in Figure (10) could not be artifacts produced by differences in 
the image quality. A radial color gradient in the observed sense has been observed for some 
active Jupiter family comets. For example, P/Tempel 2 had quite different nucleus and coma 
colors of V-R = 0.55±0.03 and 0.36±0.03, respectively (A'Hearn et al. 1989, Jewitt and Luu 
1990). Relatively blue coma is plausibly interpreted as an effect of non-geometric scattering 



in optically small particles, so perhaps Figure [10] simply shows that the particles ejected 
from the active Centaurs are, on average, small enough for geometric effects to imbue them 
with a color different from the underlying nucleus material. More complex spatial gradients 
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have been previously reported in 29P/SW1, where they may be associated with particle size 
sorting by solar radiation pressure (Cochran et al. 1982). Unfortunately, no detailed study 
of color gradients in cometary comae has been published, and the relation between the coma 
color and the color of the underlying nucleus cannot be specified with confidence. 

In summary, the active Centaurs tend to lack the ultrared matter observed in the Kuiper 
belt and some inactive Centaurs. We find no evidence that the colors of the active Centaurs 
are bimodally distributed (Figure [9]) but the number of well-measured objects is small. The 
colors of the active Centaurs might reflect a true lack of ultrared matter on their nuclei, 
or simply an artifact of color dilution with neutral-blue coma particles in the near-nucleus 



environment (Figure 10). In view of the complications imposed by coma contamination of 
the near-nucleus photometry, we are reluctant to draw any strong conclusions based on the 
correlation between Centaur color and activity. 



3.4. Orbital Distribution of the Active Centaurs 



Figure [TT] shows the distribution of the Centaurs in the semimajor axis, a, vs. eccen- 
tricity, e, plane. Active Centaurs, including those observed in the present sample plus others 
reported in the literature as being active, are shown as blue circles, while all other Centaurs 
are shown in red. Also plotted (green circles) are the Jupiter family comets (JFCs) defined 
as having a Tisserand parameter with respect to Jupiter 2 < Tj < 3. Lines are plotted in 



Figure 11 to mark the loci of orbits having perihelion distances equal to the semimajor axes 
of the giant planets. Objects above the line marked "gj", for instance, have perihelion inside 
Jupiter's orbit. 



Figure [IT] shows that the distribution of the active Centaurs in the a vs. e plane 
is different from that of the Centaurs as a whole. The active Centaurs tend to possess 
smaller than average perihelia, with most located in a band between the orbits of Jupiter 



(curve labelled qj in Figure 11) and Saturn [fis)- The median (mean) perihelion distance 
of the active Centaurs in our sample is 5.9 AU (7.6±0.9 AU, A = 9 data points) whereas 
the corresponding quantities for the inactive Centaurs in our sample are 8.7 AU (10.2±1.1 
AU, A = 11). The Centaurs as a whole have median (mean) perihelion distances 12.4 



AU (12.5±0.5 AU, A = 89). The difference is emphasized in Figure 12, which shows the 
cumulative distributions of the perihelion distance for the active and inactive Centaurs and 
for the total Centaur population. According to the Kolmogorov-Smirnov test, the statistical 
likelihood that the active and inactive Centaur perihelion distance distributions could be 
drawn by chance from the same parent population is 0.03, corresponding to about a 2.2cr 
significance. The corresponding numbers for the active Centaurs and the total Centaur 
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population are 0.01 and ~2.7cr. 

Could these differences be artifacts of an observational bias? For example, is it possible 
that Centaurs of all perihelion distances are active but that the mass-loss can be detected only 
in those with the smallest g? To investigate this possibility, we must understand the relation 
between the detection of a coma and the rate of mass-loss from the nucleus responsible for 
that coma. 



3.5. Mass-Loss Rates 

We seek to relate the optical photometry to the mass production rates of the Centaurs. 
To do this, we employ a three-step procedure. First, the photometry is used to estimate the 
scattering cross-section of the dust in the coma of each object. Second, the cross-section is 
related to the dust mass through the adoption of a size distribution of dust particles that 
is based on independent measurements from comets near 1 AU. Third, the mass-loss rate is 
inferred by dividing the dust mass by the time of residence of the dust grains as they travel 
across the projected photometry aperture. There are many unknowns in each step of this 
procedure and the derived mass-loss rates are, at best, crude estimates of the true values. 
Nevertheless, in the absence of more information, this is the best we can do. 

Magnitude (Equation [T]) is related to the geometric albedo, pr, and total cross- 
section of the coma dust particles within the projected annulus, Cd, by the inverse square 
law, conventionally expressed (Russell 1916) as 

PR<^{a)Cd = 2.25 X 10'2^i?LALlO-°-^(™'^-'"®^ (3) 

Here, Rau and Aau are the heliocentric and geocentric distances expressed in AU, a is the 
phase angle and rriQ is the apparent magnitude of the Sun viewed from 1 AU. The quantity 
$(a;) represents the phase darkening. We further take pn = 0.1, consistent with measure- 
ments of active comets in backscattering geometry (Kolokolova et al. 2004). Equations [l] and 
|3] together give the dust cross-section for each Centaur, which we list in Table |4j The dust 
cross-section and the dust mass are related through the size distribution and other param- 



eters of the dust (see Appendix). We employ Equation (A-4) to relate Cd, the cross-section 



derived from the photometry to Md, the dust mass. Values of Md are also listed in Table |4| 

The time of residence of the dust grains, at heliocentric distance Rau, in the annulus 
defined by 0i and 02 is just 
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1.5 X 10"A,„(</.2 - (j)i)/viRa 



(4) 



where Acj) = 02 — 0i, expressed in radians, and v{Rau) is the radial outflow speed of the 
dust grains from the nucleus. The heliocentric variation of the dust velocity is uncertain 
but is known to decrease with increasing heliocentric distance, albeit slowly. For example, 
measurements of expansion speeds in C/Hale-Bopp over the distance range 4 < Rau < 14 
(Biver et al. 2002) are compatible with 



viRa 



Rq 



1/4 



(5) 



where Vq = 550 m s~^ is the velocity at reference distance Rq = 5 AU. Equation ([s]) likely 
gives an upper limit to the true speed, since the Hale-Bopp measurements refer to gas. 
Entrained dust will, in general, not be perfectly coupled to the gas and therefore should 
travel more slowly. 

From Equations ^ and (|5]), the time taken for dust grains traveling radially to cross 
the annulus is written 



T{Ra 



1.5 X 10 



11 A« 



^^0 



Rn 



Rq 



1/4 



(6) 



This gives a simplistic estimate both because the grains do not radially cross the projected 
photometry annulus and because the grains presumably occupy a broad distribution of sizes, 
with the larger grains less efficiently coupled to the expanding gases of the coma and so 
traveling more slowly than small grains. Still, as a first approximation. Equation ^ provides 
a useful estimate. Column 7 of Table |4] shows that, for our data, values of r typically fall in 
the range (1 to 4)xl0^ s, or about an hour. Model mass-loss rates computed from Equations 
( |A^ and Q are listed in Table |4| 



The derived mass-loss rates are plotted against heliocentric distance in Figure [13} with 
filled circles to represent the active Centaurs and inverted triangles marking 3a limits to the 
mass-loss for those Centaurs in which no coma was detected. In passing, we note that the 
mass-loss rates for 29P estimated from our dust measurements (~5000 kg s~^) are of the 
same order as rates estimated independently, although at different times, from measurements 
of gas in this object (~1000 to 2000 kg s~^: Senay and Jewitt 1994, Crovisier et al. 1995). 
Mass loss rates from P/2004 Al have also been independently estimated from imaging data 
by Mazzotta et al. (2006). Their estimates (~100 kg s~^) are within a factor of two of 
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those in Table |4| and Figure 13 (~200 kg s~ ). We are aware of no other mass-loss rate 
estimates for the active Centaurs, but conclude from these limited comparisons that our 
simple photometric method produces values in reasonable accord with others. 



The solid lines in Figure 13 show mass-loss rates varying as 

dm 

— = 100 
dt 



Rn 



(7) 



where index n = 2, 3 and 4 and the normalization to 100 kg s"^ at Rau = 5 is arbitrary. 
The Figure shows that those Centaurs observed to be active at Rau < 13 would also be 
measurably active at larger distances, if their sublimation followed an n = 2 dependence on 
Rau- For example 166P, when displaced from 8 AU to 20 AU, would lose mass at about 60 kg 
s~^, according to the Figure, which is an order of magnitude larger than the detection limits 
estabhshed for other Centaurs (e.g. 52975, see Table ^ at this distance. [The one exception 
is the very weakly active 39P/Oterma: if displaced from 5.5 AU to larger distances an 
equilibrium coma around this object would quickly become unobservable] . Therefore, we 
conclude that the mass-loss rates in active Centaurs drop with heliocentric distance faster 
than R~^. This is significant because n = 2 is the index expected from the equilibrium 
sublimation of an exposed supervolatile ice (because essentially all the incident Solar energy 
is used for sublimation, leaving none for thermal re-radiation). Solid carbon monoxide (CO), 
for instance, would sublimate in proportion to R~^ if exposed to sunhght. This suggests 
that the absence of coma around Centaurs with Rau > 13 is not an artifact of observational 
selection acting upon a set of Centaurs losing mass in proportion to R~^. A radial variation 



with n > 3 can fit the data in Figure p!3] but this implies that the mass-loss is not driven simply 
by the sublimation of an exposed supervolatile. It is interesting to note that observations of 
comet C/Hale-Bopp as it moved out through the giant planet region of the Solar system lend 
independent support to this conclusion. The mass-loss rate from the comet closely followed 
Rau (Jewitt et al. 2008) and, indeed, coma has been consistently observed in optical images 
out as far as 26 AU (Szabo et al. 2008). 

On this basis we conclude that the activity in the Centaurs is not driven by the equi- 
librium sublimation of CO or another supervolatile. Either the mass-loss is driven by a less 
volatile material (which would have n > 2 in the observed range of heliocentric distances) or 
the restriction of the activity to smaller heliocentric distances is the result of a trigger which 
prevents mass-loss when far from the Sun. 
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4. Discussion 
4.1. Origin of the Activity 

The prime observational clues to the origin of Centaur activity are: 

» Perihelion distribution. The median perihelion distance of the active Centaurs is 
5.9 AU, whereas the median perihelion distance for the Centaurs as a whole is 12.4 



AU (Figure 12). As discussed above (see Figure 13), it is unlikely that this difference 
is the result of observational selection acting on a set of active Centaurs in the strong 
sublimation limit. Instead, the data suggest that the activity is triggered or driven by 
a temperature-related process. 

» Detection of CO. Centaur 29P is a strong source of carbon monoxide, with a sus- 
tained outgassing rate 1000 kg s~^ to 2000 kg s"^ (Senay and Jewitt 1994; Crovisier 
et al. 1995, Gunnarsson et al. 2008). CO has also been reported (but not confirmed) 
with a similar outgassing rate in Centaur (2060) Chiron (Womack and Stern 1999). 
While other Centaurs show no spectral evidence for CO sublimation or, indeed, for any 
other volatile (Bockelee-Morvan et al. 2001, Jewitt et al. 2008), these detections focus 
attention on the likely role of this supervolatile in driving mass-loss from the active 
Centaurs. 

» Line Profiles of CO. The rotational line profiles of CO in 29P are asymmetric, 
with a narrow, blue-shifted component (Senay and Jewitt 1994, Crovisier et al. 1995, 
Gunnarsson et al. 2008). This shows that the release of the CO is diurnally modulated 
and therefore that the CO source is in good thermal contact with (i.e. close to) the 
physical surface of 29P. The source cannot be deep. 

We now expand upon these points. 



Figure [M] shows solutions to the energy balance equation: 

Lo(l-A) 



47ri?2 



X 



^ dm 
eaT' + L{T) — 



in which Lq = 4x 10^^ W is the luminosity of the Sun, R is the heliocentric distance in meters, 
A is the Bond albedo, e is the emissivity and a = 5.67x10"^ W m~^ K~'^ is the Stefan- 
Boltzmann Constant. The second term on the right hand side accounts for sublimation at 
the rate dm/dt [kg m~^ s~^], and L{T) [J kg~^] is the (temperature dependent) latent heat 
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of sublimation. Conduction is neglected from Equation (|8j) based on the extremely small 
thermal inertias of two measured Centaurs (Fernandez et al. 2002) and of the nucleus of 
P/Tempel 1 (Groussin et al. 2007). Quantity x is a dimensionless parameter intended to 
represent the way in which the energy from the Sun is distributed over the surface of the 
body and is affected by both the thermal parameters of the nucleus and its rotational state. 
Physically, x may be interpreted as the ratio of the area from which absorbed sunlight is 
radiatively lost into space to the area over which energy is absorbed. We considered two 
limiting cases to bracket the likely range of x at any given heliocentric distance. A flat 
surface exposed normally to the Sun absorbs and radiates from the same area and so has x 
= 1. This limit would apply, for example, to the subsolar point of a slowly rotating body or 
to any nucleus whose projected pole direction passed through the Sun. At the other extreme, 
an isothermal sphere of radius r absorbs over vrr^ but radiates from 47rr^, and so has x = 4. 
We took L(T) from Brown and Ziegler (1979) and assumed A = 0.1, e = 0.9 (the particular 
values assumed do not have a critical effect on the computed temperatures or mass-loss rates 
provided A 1 and e 0). 



The solutions to Equation (|8]) are shown in Figure 14 for crystalline water ice (blue 
lines), CO2 ice (green lines) and CO ice (red hues). For each ice we plot two curves to show 
the low and high x values. Furthermore, at each heliocentric distance we calculated the 
critical grain size, ac, above which the gas drag force (proportional to dm/dt) is insufficient 
to launch a grain against the gravitational attraction to a nucleus of density p = 10^ kg 



m ^ and radius = 25 km. Figure 14 shows solutions of Equation M for which Oc > 1 /xm 



Grains much smaller than 1 /im are inefficient optical scatterers with characteristically blue 



colors that contradict the measured colors of the Centaurs. Figure 14 shows that CO and 
CO2 sublimation is able to eject 1 /im grains over the full range of heliocentric distances 
swept by the Centaurs while H2O sublimation can only do so for Rau < 6.5 and 8.5, in the 
low and high temperature limits, respectively. 



Figure [14] shows that none of the ices provide a convincing match to the Centaur activity. 
The Centaurs are too cold, even at their subsolar points, for water sublimation to be effective 
except close to Jupiter's orbit. For example, 2003 QD112, 166P, 167P and 2006 SX368 are 
active at 11.6, 8.6, 12.2 and 12.4 AU, all far beyond the water sublimation zone. Other 
Centaurs such as C/2001 MIO, P/2004 Al, 29P are active at distances (6.06-9.2 AU, 6.1 and 
5.8 AU, respectively) where water ice can weakly sublimate, but a comparison of the specific 
mass-loss rates with the actual mass-loss rates (Table |4]) shows that implausibly large areas 
of exposed water ice are required. More importantly, spectroscopic observations clearly show 
that the mass-loss is controlled by large fluxes of CO, at least for 29P (Senay and Jewitt 
1994, Crovisier et al. 1995). 
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Conversely, CO and CO2 ices are so volatile that they should sublimate strongly through- 
out the planetary region. If free CO ice were present on the Kuiper belt objects and Centaurs 
we should observe comae around these bodies at all heliocentric distances out to, and beyond, 
the orbit of Neptune. Since this is clearly not the case, we conclude that free CO is not re- 
sponsible for the observed activity. Non-detections of CO rotational lines have been used to 
place empirical limits to the possible fractional areas of exposed CO on Centaurs and KBOs 
near 0.1 - 1%, implying a strong depletion (Jewitt et al. 2008). Physically, this conclusion is 
easy to understand. Carbon monoxide ice sublimates so rapidly that any exposed CO could 
survive only briefly. For example, even in the cold limit, CO sublimates at 30 AU at a rate 



near dm/dt = 10~^ kg m~"^ s~^ (Figure 14), corresponding to recession of the ice surface (for 
density p = 1000 kg m^'^) at 10~^ m s~^. In one orbital period of ~150 yrs, 50 meters of 
CO would sublimate away. In 1000 orbits (only ~150kyr) the sublimation distance would 
rival the radius of the Centaur. Clearly, free CO ice cannot persist near the surface of any 
Centaur or Kuiper belt object and might only be found, if anywhere, at the bottom of deep, 
self-shadowing pits created by sublimation. The only exception to this statement applies 
to the largest Kuiper belt objects, such as Pluto, on which sublimated CO can be recycled 
through a thin, gravitationally bound atmosphere. The much smaller objects that dominate 
the Centaur population must all be free of exposed supervolatile ices on their surfaces. 

Could the CO be buried beneath a refractory layer and heated only slightly by the 
Sun? In this way, mass-loss at large Rau might be strongly suppressed relative to the n = 
2 heliocentric dependence expected of supervolatile ice. The main argument against this 
possibility is provided by the profiles of the CO rotational lines recorded in submillimeter 
spectra. The lines show a narrow, blue-shifted component that indicates that CO is released 
from the source primarily on the sun-facing side. Given that all small bodies are spinning, 
with characteristic periods typically of order tg ~ 10 hr, the CO line asymmetry requires 
that the CO be very close to the physical surface. To see this, consider a material in which 
the thermal diffusivity is k = 10~^ m^ s~^, as representative of a finely powdered dielectric. 
The diurnal thermal wave has a characteristic scale i ~ {ntsY^'^. Substituting, we estimate 
i ~ 0.06 m (6 cm). If the CO is buried any deeper than i, the heating of the CO will be 
unable to respond to the day-night insolation cycle, and the CO line profile will not show the 
sunward peak that is observed. This argument seems robust in 29P and proves that the CO 
is very close to the surface of this nucleus. However, we need more line profile measurements 
in other Centaurs in order to be sure that sunward outgassing is a general property of these 
objects. 
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4.2. Amorphous Ice in the Centaurs 

Amorphous water ice is thermodynamically unstable and converts exothermically to the 
crystalline form at a rate that depends on the ice temperature. The energy released per unit 
mass of pure water ice is about 9x10^ J kg~^, and the timescale for the conversion is given 
in years by 

(9) 

where is the activation energy, k is Boltzmann's Constant and Ej^/k = 5370 K (Schmitt 
et al. 1989). Setting Tcr = I hr (~10~'^ yr), to correspond to the timescale of a typical 
laboratory experiment, the crystallization temperature is given by Equation ^ as Tcr = 
140 K. In Solar system bodies where the timescales are much longer, significant crystalliza- 
tion can occur at correspondingly lower ice temperatures. At temperatures T < 77 K, the 
crystallization time exceeds the age of the Solar system. Amorphous ice is highly porous 
and can contain within its structure a large number of atoms and molecules trapped at the 
time of condensation of the ice from vapor. Experiments show that the trapped component 
is released upon crystallization, as the amorphous structure rearranges into a crystalline 
lattice and the trapped molecules are squeezed out (Notesco et al. 1996; 2003). On labora- 
tory timescales this happens at ~140 K but in Solar system bodies the gas release would 
occur more slowly at lower temperatures. The timescale given in Equation ^ is subject to 
uncertainties of extrapolation (the timescale is measured at temperatures much higher than 
relevant in the outer Solar system) and might be very different for ice formed in the presence 
of contaminants or particles, like silicate dust, that could act as crystallization nuclei. 

The possible role of crystallization in the activity of comets has been studied in great 
detail (e.g. Enzian et al. 1997; Klinger et al.l996; Notesco et al. 1996; 2003; Prialnik 1997). 
By providing an extra energy source, crystallization enables these models to be broadly 
successful in explaining mass-loss from comets at temperatures too low (and heliocentric 
distances too large) for crystalline water ice to sublimate. The crystallization models can be 
criticized on several levels. Disconcertingly, there is no direct evidence for the existence of 
amorphous ice in cometary nuclei (however, in fairness we note that observations provided 
no direct evidence for water ice in the comets for many decades after the formulation of 
Whipple's (1950) dirty ice nucleus model). The results of crystallization models depend 
upon a large number of unknown or poorly constrained parameters (e.g. thermal diffusivity, 
porosity, tortuosity, bulk composition, optical parameters including albedo and emissivity, 
and rotational properties) for each of which values must be assumed in order to fit the data. 
Kouchi and Sirono (2001) have further argued that the transition in real ice might not be 



T, 



cr - 3.0 X lO^^^exp 
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exothermic at all, because the sublimation of trapped volatiles liberated at crystallization 
would carry away all of the energy released. We feel that this objection has merit within a 
mean free path of the escape surface but, in buried ice, the volatiles cannot escape directly 
to space and the energy of the phase transition would remain trapped. 

While observational evidence for amorphous ice in the Centaurs is, at best, indirect, 
evidence for water ice clathrate is entirely lacking. The clathrates trap guest molecules 
within the periodic lattice structure of crystalline water ice. Hence, the volatility of clathrate 
ice is essentially the same as that of unclathrated, crystalline water ice (Klinger et al. 1996). 
The perihelion distribution of the active Centaurs, which extends far beyond the range of 



distances over which crystalline water ice can sublimate (Figure 11) is therefore inconsistent 
with the sublimation of clathrated water ice. Furthermore, the conditions of low temperature 
and low pressure likely to be found in the protoplanetary disk of the Sun, where the Centaurs 
formed, are more conducive to the formation of amorphous forms than they are to clathrates 
(Devlin 2001). If clathrates exist in the Centaurs, they have nothing to do with the observed 
outgassing activity. 

Here, we suggest a simple estimate of where in the Solar system, if exposed amorphous 
ice is present, the phase transition may be important by comparing the crystallization time, 
Tcr, to the orbital period, tk- We reason that this gives a highly conservative criterion, because 
we can be sure that amorphous ice will have crystallized on objects for which r^r <^ but 
crystallization may also have occurred on a longer timescale, up to a maximum set by the 
evolution of the orbit of the Centaur. In the spirit of the criterion, we consider circular 
orbits, for which tk = B?Ju [yr]. 



We bracket the range of likely surface temperatures as in Section (5) using Equation |8] 
with dm/dt = 0. The highest surface temperature should be that at the subsolar point, for 
which X = 1) giving T{Rau) = Ti/uli^ with Ti = 392 K. The lowest surface temperature is 
attained using the isothermal blackbody approximation, for which x = 4 and Ti = 280 K. 
Then, from Equation ^ the critical distance inside which exposed ice must be crystalline 
is found from the solution to 



llnR^^ = -A7.2Q+^^^ (10) 
Solving for 280 < Ti < 392 K gives 6.8 < Rau < 14.0 for the distances inside which the 



phase transition is expected (see also Figure 15). All of the active Centaurs have perihelia 



small enough for the crystallization of amorphous water ice to be a contributing factor. This 
observation is suggestive, albeit indirect, evidence for the existence of amorphous ice in the 
Centaurs. 



- 21 - 



There are obvious difficulties in translating the results of laboratory experiments con- 
ducted typically on timescales of hours to the relevant cometary timescales, measured in 
millions of years and longer. Still, the overall picture is very attractive. The Centaur pre- 
cursors in the Kuiper belt contain amorphous ice, laden with trapped molecules including 
CO. Once deflected into the planetary region, a thermal wave of growing amplitude begins 
to propagate beneath the surface. Regions within the Centaurs that attain peak tempera- 
tures above Tcr begin to convert to the crystalline phase, releasing trapped gases and the 
energy of the transition. This occurs preferentially on the small-perihelion Centaurs but 
post-perihelion activity can be expected because of the finite timescales for the transport of 
heat in porous media. 

Most of the active Centaurs were discovered soon before perihelion and were active at 
discovery. This is shown in Table [7| where we list the date of perihelion, the date of the 
first reported activity, the lag between these two dates, AT, the orbital period, P, and the 
ratio AT /P. The interval between these dates, AT, is typically measured in months, only 
a small fraction of the orbital period, also listed in the Table. Discovery near perihelion is 
an artifact of observational bias in magnitude limited surveys, and holds no significance for 
us here. The observation that mass-loss is underway close to and even before perihelion is 
more interesting. It tells us that there is in general no great thermal lag between the period 
of maximum heating and the on-set of mass-loss. The thermally driven process behind the 
production of coma is located close to the surface, just as we concluded earlier from the 
blue-shifted CO emission lines observed in 29P. 

Observational evidence relating to the physical state of ice on the Centaurs is limited. 
(10199) Chariklo (formerly 1997 CU26) shows water absorptions at 1.5 /im and 2.0 //m but 
there is no evidence for the 1.65 fim feature that is characteristic of crystalline ice (Dotto 
et al. 2003). Given the quality of the spectrum and the muted depths of the 1.5 and 
2.0 /im features, though, the absence of the 1.65 /im band cannot be used as evidence for 
amorphous ice. (2060) 95P/Chiron likewise shows water ice absorptions (Foster et al. 1999, 
Luu et al. 2000) but with absorption band strengths that are modulated by the activity of 
the coma. The bands are again too weak and the data are too poor to permit us to rule 
definitively on the crystalline vs. amorphous state of the ice. The activity in 2060 Chiron 
has been modeled as a product of crystallization heating by Prialnik et al. (1995). 

4.3. Amorphous Ice in the Comets and Kuiper Belt Objects 

The available data are consistent with the presence of amorphous ice in the Centaurs. 
What does this mean for the dynamically related populations of Jupiter family comets ( JFCs) 
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and Kuiper belt objects? 

A fraction of the Centaurs will eventually be captured as JFCs, suggesting that these 
bodies should also contain amorphous ice. Observationally, direct measurements of the ice in 
comets are very few, and the crystalline state is not well measured. However, as noted above, 
crystallization has long been invoked by modelers as an important process in the JFCs, with 
an apparent success that is tempered by knowledge of the large number of poorly constrained 
parameters that must be assumed. Crystallization models also provide plausible fits to the 
activity of some long-period and dynamically new (Oort cloud) comets (Meech et al. 2009). 
A fair statement would be that there is no direct evidence concerning the crystallinity of 
cometary water ice and nothing to strongly contradict the presence of ice in the amorphous 
state. 

A somewhat different circumstance prevails for the precursors to the Centaurs located in 
the Kuiper belt. There, near infrared spectral observations provide a seemingly contradictory 
result; the 1.65 //m absorption band which is diagnostic of the presence of crystalline water 
ice is clearly observed in many KBOs (Jewitt and Luu 2004, Takato et al. 2006, Trujillo et 
al. 2007). Some fraction of the ice could be amorphous and still fit the measured profile 
of the 1.65 /im band (Merlin et al. 2007, Pinilla-Alonso et al. 2009). On the other hand, 
there are no well-observed objects in which the characteristic 1.5 yum and 2.0 /im bands of 
water ice are observed but the 1.65 \xva. crystalline ice band is absent. Therefore, where the 
data are good enough to make a diagnosis, the available spectra provide no evidence that 
amorphous is ubiquitous on the KBOs. If the Centaurs contain amorphous ice, and if they 
recently escaped from long-term storage in the Kuiper belt, it is logical to conclude that 
ice in the KBOs should be, at least in part, amorphous. Certainly, amorphous ice could 
survive indefinitely in the KBOs given the very low radiation equilibrium temperatures in 
the Kuiper belt. How might this seeming contradiction be explained? 

One possibility is that this is a size-selection effect. The KBOs bright enough for the di- 
agnostic 1.65 yun band to be spectroscopically detected tend to be large (~1000 km) objects. 
Heating from their gravitational binding energy (McKinnon 2002) or from trapped radioac- 
tive nuclei may have caused the ice in these large bodies to have crystallized, so explaining 
the spectroscopic data (Merk and Prialnik 2003). Smaller bodies of size comparable to the 
Centaurs listed in Table [s] (radii of a few x 10 km or less) could have escaped crystallization 
but are effectively too faint to observe spectroscopically in the infrared. A second possibil- 
ity is that the KBOs and Centaurs are stratified, with crystalline ice at the visible surface 
and amorphous ice beneath. For example, some process such as heating by micrometeorite 
bombardment could crystallize surface ice while leaving deeper ice in the amorphous state. 
Evidently, these two possibilities are not mutually exclusive. 
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4.4. The Future 

A number of exciting observational prospects exist for making progress in this subject. 
First, the assertion that the active Centaurs are powered by the crystalhzation of amorphous 
ice is open to an observational test: this mechanism could be refuted if future observations 
reveal outgassing activity in Centaurs having perihelion distances much larger than the crit- 



ical range for crystallization, as given by Equation (10). For example, the discovery of active 
Centaurs with perihelion q > 20 AU would be difficult to reconcile with the crystallization 
of amorphous ice, unless the dynamics showed that they had been very recently scattered 
outwards from an orbit with a smaller perihelion. Second, deep rotational line spectra with 
ALMA should show whether CO outgassing drives the mass-loss in other active Centaurs, 
as is expected in the present hypothesis. Rotational line spectra having high sensitivity will 
also allow us to measure the mass-loss rates with greater accuracy, instead of relying on im- 
perfect estimates based on observations of the dust, as in this work. Third, measurements of 
the CO line profiles will show whether or not the mass-loss is mainly sunward (and therefore 
close to the physical surface), as is observed in 29P/Schwassmann-Wachmann 1. 

Future observations might also show why some Centaurs remain inactive even when their 



perihelia are small. In Figure (11), more than half of the Centaurs with perihelia between 



the orbits of Jupiter and Saturn appear inactive, although all are, presumably, hot enough 



to satisfy the crystallization criterion in Equation (10). This could be because the mass- 
loss is intermittent or "bursty", on timescales longer than the orbital period, as a result of 
instabilities in the crystallization front. Alternatively, these hot but inactive Centaurs might 
have been devolatilized in their near-surface regions as a result of past activity (i.e. they 
could have been in their current orbits for longer, as a group, than the active Centaurs). 
Still another possibility is that some Centaurs are deficient in volatiles, or perhaps lacking 
in the amorphous ice whose exothermic transformation appears to provide a plausible gas 
and energy source for the active Centaurs. Physical studies, careful long-term photometric 
monitoring and dynamical investigations of the Centaurs will help us to understand which 
of these possibilities might be correct. 
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5. Summary 

We present new observations of the active Centaurs, defined as Solar system bodies 
whose orbits have both periheUa and semimajor axes between the orbits of Jupiter (at 5.2 
AU) and Neptune (30 AU) and which are not in 1:1 resonance with the giant planets. 

1. Including observations from the present work and from the literature, about 13% (12 of 
92) Centaurs are measurably active, as shown by the presence of dust comae in optical 
images. 

2. The active Centaurs have statistically smaller perihelion distances (median 5.9 AU) 
than the inactive Centaurs in our sample (median 8.7 AU) and smaller than the me- 
dian perihelion distance of all known Centaurs (12.4 AU). Simulations show that the 
difference is unhkely to be an artifact of observational selection. Instead, the smaller 
distances of the active Centaurs appear to reflect the action of a thermally-driven 
trigger or activity source. 

3. The surface temperatures of the active Centaurs are too low (and the inferred mass-loss 
rates are too large) for the sublimation of crystalline water ice to constitute a plausible 
source of the activity. 

4. Conversely, although gaseous carbon monoxide (CO) has been reported in two Centaurs 
(29P and Chiron), CO ice is too volatile to explain the measured perihelion distribution. 
If present, near-surface CO ice would strongly sublimate all the way out to Neptune's 
orbit, whereas no active Centaurs are known so far from the Sun. 

5. The active Centaurs have orbital periods which are long compared to the timescale 
for the crystallization of water ice at their distance. No active Centaurs in which this 
inequality is reversed have yet been found. This fact is consistent with (but does not 
prove) cystallization as the trigger, or driver of activity. 

6. If Centaurs contain amorphous ice, as suggested by the available data, then their 
precursor bodies in the Kuiper belt must likewise be at least partly amorphous and 
their progeny, the nuclei of Jupiter family comets, should also contain amorphous ice. 

7. Most active Centaurs are slightly redder than the Sun at optical wavelengths, but they 
appear devoid of the ultrared matter found in the inactive Centaurs and Kuiper belt 
objects. Physical interpretation of the colors, however, is stymied by coma contamina- 
tion in the active Centaurs. 
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APPENDIX 



The connection between the cross-section of the particles and their mass depends largely 
upon the size distribution of the particles, but also on their density and scattering efficiency. 
For spherical particles of a single radius, a, the mass, M^, and the cross-section, Cd = vra^, 
are (in the geometric optics limit) proportional: 



Md = -paCd (A-1) 

where p is the grain bulk density. 

This simple relation may be generalized for power-law size distributions in which the 
number of particles with radius between a and a + da is Ta~'^da, where F and q are constants 
of the distribution and the relation holds between minimum and maximum particle radii a_ 
and a+, respectively. The combined cross-section of such a distribution is 



Cd = T^T I a^-'^da (A-2) 

while the total mass is 



Md 



-vrF 



1^ '^da 



(A-3) 



Results obtained from the study of active comets show that the dust does not exactly 
follow a power law size distribution, but that a reasonable approximation is obtained by 
setting q = 3.5, a_ = 0.1 pm and a+ = 1 cm, (e.g. Griin et al. 2001). Comparing Equations 



(A-2) and (A-3) then shows that the mass is given by 
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M, = ^p{a_a+f'C, (A-4) 

and this is the relation we employ here. The size of the average scatterer in this case is a 
= (a_a+)^/^ ~ 30/im, consistent with the observation that the comae of the Centaurs are 
neutral or reddish compared to the Sun, not blue (as would be expected if a -C A, where 
A ~ 0.5 fim is the wavelength of observation). 



Other relations can be trivially computed from Equations (A-2) and (A-3) and still 
more involved estimates can be obtained by including the wavelength- and size-dependent 
scattering efficiencies of the particles. However, given our basic ignorance of the properties 
of the dust in Centaurs, these extra steps hardly seem appropriate. In fact, uncertainties in 
the dust parameters set a fundamental limit to our ability to estimate dust mass production 
rates from optical data and the absolute values of the mass-loss rates derived from Equation 



(A-4) could be wrong by an order of magnitude or more. In this paper, we use Equation 
(A-4), with (a_a+)^/^ = 30/im and p = 1000 kg m~^, simply as a useful metric with which 



to compare relative mass-loss rates amongst the Centaurs. 
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Table 1. Journal of Observations 





TTT rin+fi 
U ± UdiLe 


±ei 


D a 


A b 


[deg]*' 


C/2001 MIO 


2002 Sep 08 


UH 2.2-m 


6.056 


5.349 


4.16 


C/2001 MIO 


2004 Oct 10 


Keck 10-m 


9.147 


8.291 


3.39 


P/2004 Al(LONEOS) 


2006 Aug 01 


UH 2.2-m 


6.079 


5.961 


9.62 


39P/Otcrma 


2002 Sep 07 


UH 2.2-m 


5.484 


4.684 


6.95 


29P/SW1 


2002 Sep 08 


UH 2.2-m 


5.810 


4.991 


6.24 


174P/Echeclus (60558) 


2006 Feb 25 


Keck 10-m 


12.968 


12.230 


3.00 


P/2005 T3 (Read) 


2007 Sep 19 


UH 2.2-m 


6.442 


6.368 


8.97 


P/2005 T3 (Read) 


2007 Oct 12 


Keck 10-m 


6.459 


6.026 


8.27 


P/2005 S2 (Skiff) 


2007 Oct 12 


Keck 10-m 


6.552 


5.866 


6.70 


165P/Lincar 


2003 Jan 08 


Keck 10-m 


8.300 


8.123 


6.75 


(63252) 2001 BL41 


2002 Feb 22 


UH 2.2-m 


8.456 


7.547 


2.79 


(63252) 2001 BL41 


2002 Feb 23 


UH 2.2-m 


8.458 


7.555 


2.89 


2000 GM137 


2001 Feb 16 


UH 2.2-m 


6.970 


6.348 


6.62 


2000 GM137 


2001 Feb 17 


UH 2.2-m 


6.970 


6.336 


6.54 


2003 QD112 


2004 Oct 10 


Keck 10-m 


11.625 


10.692 


1.82 


(145486) 2005 UJ438 


2007 Oct 12 


Keck 10-m 


9.280 


8.609 


4.73 


(32532) Thereus 


2007 Feb 20 


Keck 10-m 


10.792 


11.233 


4.60 


166P/2001 T4 


2001 Nov 18 


UH 2.2-m 


8.590 


7.841 


4.49 


166P/2001 T4 


2002 Sep 07 


UH 2.2-m 


8.573 


7.937 


5.43 


(5145) Pholus 


2004 Feb 17 


Keck 10-m 


18.538 


18.488 


3.05 


(148975) 2001 XA255 


2002 Dec 07 


Keck 10-m 


15.062 


14.262 


2.24 


(148975) 2001 XA255 


2003 Jan 09 


Keck 10-ni 


14.965 


13.982 


0.14 


167P/2004 PY42 


2004 Oct 10 


Keck 10-ni 


12.234 


11.840 


4.36 


2006 SX368 


2007 Sep 19 


UH 2.2-m 


12.398 


11.703 


3.44 


(54598) Bienor 


2006 Aug 23 


UH 2.2-m 


18.369 


17.470 


1.49 


2002 PN34 


2006 Aug 23 


UH 2.2-m 


14.661 


13.678 


0.95 


2002 VG131 


2002 Dec 07 


Keck 10-m 


14.870 


14.431 


3.24 


2001 XZ255 


2003 Jan 29 


Keck 10-m 


16.079 


15.096 


0.09 


2001 XZ255 


2004 Oct 10 


Keck 10-m 


16.167 


16.385 


3.43 


(52975) Cyllarus 


2004 Oct 10 


Keck 10-m 


21.089 


20.542 


2.30 



^Heliocentric distance in AU 

b Geocentric distance in AU 
Phase angle in degrees 
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Table 2. Orbital Parameters of the Observed Centaurs^ 



Object 


q [AU] b 


a" 




ed 


i [deg]" 


Perihelion 


C/2001 MIO 


5.303 


26.660 


0, 


,801 


28.0 


2001 Jun 


21 


P/2004 Al(LONEOS) 


5.463 


7.895 


0, 


,308 


10.6 


2004 Aug 


25 


39P/Otcmia 


5.471 


7.256 


0, 


,246 


1.9 


2002 Dec 


21 


29P/SW1 


5.722 


5.986 


0, 


,044 


9.4 


2004 Jun 


30 


174P/Echeclus (60558) 


5.849 


10.740 


0, 


,455 


4.3 


2015 May 


06 


P/2005 T3 (Read) 


6.202 


7.507 


0, 


,174 


6.3 


2006 Jan 


13 


P/2005 S2 (Skiff) 


6.398 


7.964 


0, 


,197 


3.1 


2006 Jun 


30 


165P/Lincar 


6.830 


18.03 


0, 


,621 


15.9 


2000 Jun 


15 


(63252) 2001 BL41 


6.880 


9.816 


0, 


,299 


12.4 


1998 Feb 


09 


2000 GM137 


6.951 


7.907 


0, 


,121 


15.8 


2002 Feb 


24 


2003 QD112 


7.935 


18.974 


0, 


,582 


14.5 


1998 May 


29 


145486 


8.255 


17.527 


0, 


,529 


3.8 


2010 Jul 


01 


(32532) Thcrcus 


8.524 


10.615 


0, 


,197 


20.4 


1999 Feb 


15 


166P/2001 T4 


8.564 


13.880 


0, 


,383 


15.4 


2002 May 


20 


(5145) Pholus 


8.685 


20.292 


0, 


,572 


24.7 


1991 Sep 


24 


(148975) 2001 XA255 


9.387 


30.183 


0, 


,689 


12.7 


2010 Jun 


04 


167P/2004 PY42 


11.784 


16.140 


0. 


,270 


19.1 


2001 Apr 


10 


2006 SX368 


11.963 


22.236 


0. 


,462 


36.3 


2010 May 


16 


(54598) Bienor 


13.155 


16.485 


0, 


,202 


20.8 


2027 Dec 


19 


2002 PN34 


13.328 


30.781 


0, 


,567 


16.7 


2001 Nov 


11 


2002 VG131 


14.869 


17.492 


0, 


,150 


3.2 


2002 Nov 


09 


2001 XZ255 


15.353 


15.910 


0, 


,035 


2.6 


1984 Jun 


22 


(52975) Cyllarus 


16.286 


26.438 


0, 


,384 


12.6 


1989 Oct 


27 



^Objects are listed in the order of increasing perihelion distance 

''Perihelion distance in AU 

'^Semimajor axis in AU 

'^Orbital eccentricity 
° Orbital inclination 
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Table 3. R-Band Photometry for Coma Search 



Object 


UT Date'' 










mH(l,l,0)<= 1 


"e [km]' 


C/2001 MIO 


2002 


Sep 


08 


18. 


,85 


18 


,53 


20, 


.02±0.04 


>10 


.92 


<12 


P/2004 Al 


2006 


Aug 


01 


19. 


.76 


19, 


,44 


20, 


,95±0.04 


>11 


,33 


<10 


39P/Oterma 


2002 


Sep 


07 


21. 


,83 


21 


,75 


24, 


,64±0.05 


>14, 


.27 


<3 


29P/SW1 


2002 


Sep 


08 


16. 


,32 


15, 


,86 


17, 


,01±0.05 


>8, 


.52 


<35 


174P/Echeclus (60558) 


2006 


Feb 


25 


20. 


,01 


19, 


,69 


21, 


,19±0.05 


>8 


,70 


<33 


P/2U05 bZ (bkiit) 


2007 


Oct 


12 


22. 


,97 


23, 


,07 




>26.18 


14, 


,54 


2 


(63252) 2001 BL41 


2002 


Feb 


22 


20. 


,35 


20, 


,33 




>23.85 


11, 


.03 


11 


(63252) 2001 BL41 


2002 


Feb 


23 


20. 


,39 


20, 


,37 




>25.17 


11, 


.06 


11 


2000 GM137 


2001 


Feb 


16 


22. 


,61 


22, 


,63 




>25.62 


13, 


,88 


3 


2000 GM137 


2001 


Feb 


17 


22. 


,61 


22, 


,62 




>27.42 


13, 


,89 


3 


2003 QD112 


2004 


Oct 


10 


22. 


,32 


22, 


,01 


23, 


,55±0.04 


>11, 


,62 


<9 


(145486) 2005 UJ438 


2007 


Oct 


12 


20. 


,17 


20, 


,15 




>23.87 


10, 


,25 


16 


(32532) Thereus 


2007 


Feb 


20 


19. 


,87 


19, 


,87 




>23.39 


9, 


,05 


28 


166P/2001 T4 


2001 


Nov 


18 


20. 


,01 


19, 


,78 


21, 


,60±0.05 


>10, 


,47 


<14 


166P/2001 T4 


2002 


Sep 


07 


19. 


,65 


19, 


,27 


20, 


,63±0.05 


>10, 


.04 


<18 


(5145) Pholus 


2004 


Feb 


17 


19. 


,96 


19, 


,94 




>25.51 


7, 


.00 


72 


(148975) 2001 XA255 


2002 


Dec 


07 


22. 


,74 


22, 


,74 




>26.93 


10, 


.82 


12 


(148975) 2001 XA255 


2003 


Jan 


09 


22. 


,39 


22, 


,30 




>26.20 


10, 


.74 


13 


167P/2004 PY42 


2004 


Oct 


10 


20. 


,69 


20, 


,64 


24, 


,04±0.05 


>9 


.50 


<23 


2006 SX368 


2007 


Sep 


19 


20. 


,15 


20, 


,08 


23.8±0.1 


9, 


.00 


28 


(54598) Bienor 


2006 


Aug 


23 


21. 


,11 


21, 


,53 




>25.90 


8, 


.38 


38 


2002 PN34 


2006 


Aug 


23 


19. 


,92 


19, 


,89 




>26.23 


8, 


.26 


40 


2002 VG131 


2002 


Dec 


07 


23. 


,08 


23, 


,25 




>26.13 


11. 


.10 


11 


2001 XZ255 


2003 


Jan 


09 


22. 


,44 


22, 


,34 




>26.14 


10, 


.48 


14 


2001 XZ255 


2004 


Oct 


10 


21. 


,11 


23, 


,11 




>26.26 


8, 


,66 


34 


(52975) Cyllarus 


2004 


Oct 


10 


21. 


,67 


21, 


,66 




>26.29 


8, 


,23 


41 



^UT Date of the oberservation 

''Magnitude within a 2.2 arcsec radius aperture 

"^Magnitude within a 3.3 arcsec radius aperture 

"^Coma magnitude within the 2.2 to 3.3 arcsec annulus (Equation^ 

''Red magnitude of the nucleus, computed from 7712.2 corrected to unit heliocentric and geocentric 
distance and zero phase angle (Equation |2]| 

'Effective radius of the nucleus in km assuming red geometric albedo pn = 0.1. Limits are given 
where coma was detected. 
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Table 4. Model Mass Loss Rates 



Object 


R [AUf 


A [AU]'^ 


a [dcg]" 






r [s]t 


Md [kg]g 


dMd h 
dt 


C/2001 MIO 


6.056 


5.349 


4.16 


20.0 


5.77x10^ 


8.16x10^ 


2.31xl0'5 


2.8x10^ 


P/2004 Al 


6.079 


5.961 


9.62 


21.0 


4.55x10^ 


9.10xl03 


1.82x10*5 


2.0x10^ 


39P/Oterma 


5.484 


4.684 


6.95 


24.6 


7.61x10^ 


6.97x103 


3.04x10"' 


4.3x10° 


29P/SW1 


5.810 


4.991 


6.24 


17.0 


9.61x10** 


7.53x10^ 


3.84x10'^ 


5.1x10^ 


174P 


12.970 


12.230 


3.00 


21.2 


1.85x10** 


2.26x10"' 


7.41x10'' 


3.2x10^ 


P/2005 82 


6.552 


5.866 


6.70 


>26.28 


<2.93xl05 


9.12x10^ 


<1. 18x10"' 


<1.3xl0° 


63252 


8.458 


7.555 


2.89 


>25.17 


<1.23xl0'5 


1.25x10"' 


<4.93xl0"' 


<3.9xl0° 


2000 GM137 


6.970 


6.336 


6.54 


>27.42 


<1.24xl05 


1.00x10" 


<4.95xl03 


<0.5xl0'' 


2003 QD112 


11.625 


10.692 


1.82 


23.5 


1.33x10^ 


1.92x10"' 


5.31x10^ 


2.8x10^ 


145486 


9.280 


8.609 


4.73 


>23.87 


<6.63xl0'5 


1.46x10"' 


<2.65xl05 


<1.8xl0i 


32532 


10.792 


11.233 


4.60 


>23.39 


<1. 80x10^ 


1.98x10"' 


<7.19xl05 


<3.6xl0i 


166P/2001 T4 


8.573 


7.937 


5.43 


20.6 


1.10x10* 


1.32x10"' 


4.40x10'' 


3.3x10^ 


(5145) Pholus 


18.538 


18.488 


3.05 


>25.51 


<1. 08x10^ 


3.73x10"' 


<4.30xl05 


<1.2xl0i 


148975 


15.062 


14.262 


2.24 


>26.93 


<1. 37x10*^ 


2.73x10"' 


<5.49xl0"' 


<2.0xl0° 


167P/2004 PY42 


12.234 


11.840 


4.36 


24.0 


1.31x10^ 


2.15x10"' 


5.24x10^ 


2.4x10^ 


2006 SX368 


12.398 


11.703 


3.44 


23.8 


1.52x10^ 


2.13x10"' 


6.10x105 


2.9x10^ 


54598 


18.369 


17.470 


1.49 


>25.90 


<6. 03x10*^ 


3.52x10"' 


<2.41xl05 


<6.9xl0'5 


2002 PN34 


14.661 


13.678 


0.95 


>26.23 


<2. 11x10*^ 


2.60x10"' 


<8.45xl0"' 


<3.2xl0° 


2002 VG131 


14.870 


14.431 


3.24 


>26.13 


<3.10xl0'5 


2.75x10"' 


<1.24xl05 


<4.5xl0'' 


2001 XZ255 


16.167 


16.385 


3.43 


>26.26 


<3. 76x10*5 


3.19x10"' 


<1.50xlO"5 


<4.7xl0° 


52975 


21.089 


20.542 


2.30 


>26.29 


<7.03xl0'5 


4.28x10"' 


<2.81xl05 


<6.6xl0'' 



Heliocentric distance 
''Geocentric distance 
'^Phase angle 

■^Coma annulus magnitude from Table (jsjl. We have taken the most sensitive limit for each object in which a coma 
was not detected 

"Derived dust cross-section (Equation [sjl 

^Photometry annulus crossing time (Equation [6| 

^Derived coma mass [kg], from Cd and Equation |A-4| 

'^Derived mass-loss rate [kg s^^^], from columns (g) and (f) 
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Table 5. Color and Multi- Aperture Photometry 



Object 


UT 2005 














V-R'^ 




R-P 




29P/SW1 


2002 


Sep 08 


1, 


.1 


16, 


.56±0. 


.01 


0.79±0.03 


0, 


.51±0, 


.03 


0. 


.52±0. 


,03 


29P/SW1 


2002 


Sep 08 


2, 


.2 


16. 


.32±0. 


.01 


0.78±0.03 


0, 


.50±0, 


.03 


0. 


.52±0. 


.03 


29P/SW1 


2002 


Sep 08 


3, 


.3 


15. 


.86±0. 


.01 


0.79±0.04 


0, 


.47±0, 


.03 


0. 


.50±0. 


.03 


29P/SW1 


2002 


Sep 08 


6, 


.6 


15. 


,83±0. 


.01 


0.75±0.03 


0, 


,51±0, 


.03 


0. 


.48±0. 


.03 


39P/Otcrma 


2002 


Sep 07 


1, 


.1 


22. 


,08±0. 


,05 


0.89±0.07 


0, 


,41±0, 


.05 


0. 


,47±0, 


.10 


39P/Otcrma 


2002 


Sep 07 


2, 


.2 


21. 


,83±0. 


,05 


0.74±0.07 


0, 


,45±0, 


.05 


0. 


,59±0, 


.10 


39P/Oterma 


2002 


Sep 07 


3, 


.3 


21. 


,75±0. 


,05 


0.72±0.07 


0, 


,40±0, 


.05 


0. 


,44±0, 


.10 


166P/2001 T4 


2001 


Nov 18 


1, 


.1 


20. 


.47±0. 


,03 


0.87±0.04 


0, 


,69±0. 


.03 


0. 


,78±0, 


.04 


166P/2001 T4 


2001 


Nov 18 


2, 


.2 


20. 


.01±0. 


,03 


0.87±0.04 


0, 


.73±0. 


.03 


0. 


,71±0, 


.04 


166P/2001 T4 


2001 


Nov 18 


3, 


.3 


19. 


,78±0. 


,03 


0.87±0.04 


0. 


.66±0, 


.04 


0. 


,73±0, 


.04 


166P/2001 T4 


2001 


Nov 18 


6, 


.6 


19. 


,35±0. 


,02 


0.70±0.04 


0, 


.73±0, 


.03 




_ 




166P/2001 T4 


2002 


Sep 07 


1, 


.1 


19. 


,91±0. 


,03 


_ 


0, 


.70±0, 


.03 


0. 


,69±0, 


.04 


166P/2001 T4 


2002 


Sep 07 


2, 


.2 


19. 


,65±0. 


,03 


_ 


0. 


.63±0, 


.03 


0. 


,70±0, 


.04 


166P/2001 T4 


2002 


Sep 07 


3, 


.3 


19. 


.27±0. 


,03 


_ 


0. 


.52±0, 


.03 


0. 


,76±0, 


.04 


166P/2001 T4 


2002 


Sep 07 


3, 


.3 


19. 


.03±0. 


,02 


_ 


0. 


.52±0. 


.03 


0. 


,76±0, 


.04 


167P/2004 PY42 


2004 


Oct 10 


1, 


.1 


20. 


.80±0. 


,03 


0.77±0.03 


0. 


.55±0, 


.03 


0. 


,54±0, 


.02 


167P/2004 PY42 


2004 


Oct 10 


2 


.2 


20. 


,69±0. 


,02 


0.80±0.03 


0, 


,51it0 


.03 


0. 


,51±0, 


.03 


167P/2004 PY42 


2004 


Oct 10 


3, 


.3 


20. 


,64±0. 


,04 


0.76±0.03 


0. 


,49±0, 


.04 


0. 


,51±0, 


.04 


C/2001 MIO 


2002 


Sep 08 


1, 


.1 


19. 


,14±0. 


,02 


1.22±0.08<^ 














C/2001 MIO 


2002 


Sep 08 


2, 


.2 


18. 


.85±0. 


,02 


1.19±0.08<* 














C/2001 MIO 


2002 


Sep 08 


3, 


.3 


18. 


.53±0. 


,03 


1.17±0.08'^ 














C/2001 MIO 


2002 


Sep 08 


6, 


.6 


18. 


.52±0. 


,04 


1.16±0.08'^ 














C/2001 MIO 


2004 


Oct 10 


1, 


.1 


22. 


,76±0. 


,06 


0.77±0.06 


0. 


.35±0, 


.06 


0. 


,52±0, 


.08 


C/2001 MIO 


2004 


Oct 10 


2, 


.2 


22. 


,08±0. 


,06 


0.66±0.10 


0. 


.49±0. 


.10 


0. 


,42±0, 


.08 


174P/Echcclus (60558) 


2006 


02 25 


1, 


.1 


20. 


,27±0. 


,02 
















174P/Echeclus (60558) 


2006 


02 25 


2, 


.2 


20. 


,01±0. 


,02 
















174P/Echeclus (60558) 


2006 


02 25 


3, 


.3 


19. 


.69±0. 


,02 
















174P/Echeclus (60558) 


2006 


02 25 


6, 


.6 


18. 


.87±0. 


,02 
















2003 QD112 


2004 


10 10 


1, 


.1 


22. 


,49±0. 


,03 
















2003 QD112 


2004 


10 10 


2, 


.2 


22. 


,32±0. 


,03 
















2003 QD112 


2004 


10 10 


3, 


.3 


22. 


,18±0. 


,03 
















2003 QD112 


2004 


10 10 


6, 


.6 


22. 


■OlitO. 


,05 
















P/2004 Al 


2006 


08 01 


1, 


.1 


20. 


.35±0. 


,02 
















P/2004 Al 


2006 


08 01 


2, 


.2 


19. 


.76±0. 


,02 
















P/2004 Al 


2006 


08 01 


3, 


.3 


19. 


.44±0. 


,02 




0. 


.32±0. 


.03 








P/2004 Al 


2006 


08 01 


6, 


.6 


19. 


OlitO. 


,02 




0. 


.39±0. 


.03 









^Radius of the photometry aperture, in arcsec. 
''Apparent R-band magnitude and its la uncertainty. 
Color indices within each aperture 
''The B-R color is reported on this object only 
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Table 6. Adopted Colors of the Active Centaurs 



Object 


UT 2005 




Rb 




V-R'^ 


R-P 


29P/SW1 (coma) 


2002 


Sep 


08 


2.2 


16.32±0.02 


0.78±0.03 


0.50±0.03 


0.52±0.03 


39P/Oterma (coma) 


2002 


Sep 


07 


2.2 


21.83±0.05 


0.74±0.07 


0.45±0.05 


0.59±0.10 


166P/2001 T4 (coma) 


2001 


Nov 


18 


2.2 


20.01±0.02 


0.87±0.03 


0.73±0.03 


0.71±0.03 


(2060) 95P/Chiron<l 


2001 


Jun 


13 


1.4 


15.77±0.04 


0.60±0.03 


0.32±0.03 


0.64±0.03 


174P/Echeclus (60558) " 


2006 


Mar 


23 


1.0 


20.13±0.1 


0.71±0.15 


0.56±0.15 




174P/Echeclus (60558) ' 


2006 


Feb 


24 


30 


15.18±0.04 


0.81±0.07 


0.50±0.06 


0.58±0.05 


C/2001 MIO (coma) 


2004 


Oct 


10 


2.2 


22.08±0.06 


0.66±0.10 


0.51±0.03 


0.51±0.03 


167P/2004 PY42 (nucleus) 


2004 


Oct 


10 


2.2 


20.69±0.02 


0.80±0.03 


0.49±0.10 


0.42±0.08 


P/2004 Al 


2006 


Oct 


01 


3.3 


19.76±0.02 




0.32±0.03 




Solar Colors 












0.67 


0.36 


0.35 



''Radius of the circular aperture within which the measurement was taken, in arcseconds 
^ij-band magnitude 
Color index, as listed, with Icr uncertainty 
tJRomon-Martin et al. 2003 

''Rousselot 2008; this small aperture measurement samples the primary nucleus 

'Bauer ct al. (2008); this large aperture measurement samples the primary and secondary nuclei plus 

the associated dust envelope 
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Table 7. Onset of Activity 



Object 


Perihelion Date* 


Activity Date*" 


AT [yr]<^ 


p [yr]d 


AT 
P 


Comment 


C/2001 MIO 


2001 Jun 21 


2001 Jun 20 


0.0 


137.7 


0.00 


Active at discovery 


P/2004A1(LONEOS) 


2004 Aug 25 


2004 Jan 13 


-0.6 


22.2 


-0.03 


Active at discovery 


39P/Oterma 


2002 Dec 21 


2002 Sep 07 


-0.3 


19.5 


-0.02 




29P/SW1 


2004 Jun 30 






14.6 




Always active 


174P/Echeclus (60558) 


2015 May 06 


2006 Jan 


-9.4 


35.2 


-0.27 




P/2005T3 (Read) 


2006 Jan 13 


2005 Oct 


-0.3 


20.6 


-0.01 


Active at discovery 


P/2005 S2 (Skiflf) 


2006 Jun 30 


2005 Oct 


-0.8 


22.5 


-0.04 


Active at discovery 


165P/Lincar 


2000 Jun 15 


2000 Jan 29 


-0.5 


76.6 


-0.01 


Active at discovery 


2003 QD112 


1998 May 29 


2004 Oct 10 


+4.4 


82.7 


+0.05 




(2060) 95P/Chiron 


1996 Feb 14 


1989 Apr 10 


-6.8 


50.8 


-0.13 


Active since 1989 


166P/2001 T4 


2002 May 20 


2001 Oct 15 


-0.6 


51.9 


-0.01 


Active since discovery 


167P/2004 PY42 


2001 Apr 10 


2004 Aug 11 


+3.3 


64.8 


+0.05 





''Date of closest perihelion 

''Dates of first reported cometary activity 

'^Time lag, in years, between perihelion and first date of reported activity. Negative means activity was detected 
before perihelion 

''Orbital period, in years 



- 38 - 





29P 



Fig. 1. — 29P/Schwassmann-Wachmann 1 imaged in the R-band filter on UT 2002 Sep 08. 
This is a 300 sec integration with the University of Hawaii 2.2-m telescope. The telescope 
was tracked non-sidereally to follow the motion of the Centaur. Field shown has North to 
the bottom, East to the right, and is 130 arcsec wide. The centrally located Centaur is 
marked with a bar. 
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Fig. 2. — 39P/Oterma imaged in the R-band filter on UT 2002 Sep 08. This is the sum 
of two consecutive 900 sec integrations with the University of Hawaii 2.2-m telescope. The 
telescope was tracked non-sidereally to follow the motion of the Centaur. Field shown has 
North to the bottom, East to the right, and is 130 arcsec wide. 
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Fig. 3.— 166P/2001 T4 imaged in the R-band filter on UT 2002 Sep 07. This is a 300 sec 
integration with the University of Hawaii 2.2-m telescope. The telescope was tracked non- 
sidereally to follow the motion of the Centaur. The diagonal hne is from a passing satellite. 
Field shown has North to the bottom, East to the right, and is 130 arcsec wide. 
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Fig. 4.— 167P/2004 PY42 imaged in the R-band filter on UT 2002 Sep 08. This is the 
sum of two 200 sec integrations with the Keck 10-m telescope. The telescope was tracked 
non-sidereally to follow the motion of the Centaur. Field shown has North to the bottom, 
East to the right, and is 130 arcsec wide. 
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Fig. 5.— C/2001 MIO imaged in the R-band filter on UT 2002 Sep 08. This is a 300 sec 
integration with the University of Hawaii 2.2-m telescope. The telescope was tracked non- 
sidereally to follow the motion of the Centaur. Field shown has North to the bottom, East 
to the right, and is 130 arcsec wide. 
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Fig. 6.— P/2004 Al (Loneos) imaged in the R-band filter on UT 2006 Jul 01. This is a 300 
sec integration with the University of Hawaii 2.2-m telescope. The telescope was tracked 
non-sidereally to follow the motion of the Centaur. Field shown has North to the bottom, 
East to the right, and is 130 arcsec wide. 
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Fig. 7.— 2003 QD112 imaged in the R-band filter on UT 2004 Oct 10. This is a 300 sec 
integration with the Keck 10-m telescope. The telescope was tracked non-si dereally to follow 
the motion of the Centaur. Field shown has North to the bottom, East to the right, and is 
130 arcsec wide. The centrally located Centaur is marked with a bar. 
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Fig. 8.— 174P/Echeclus (60558) imaged in the R-band filter on UT 2005 Feb 25. This is a 
300 sec integration with the Keck 10-m telescope. The telescope was tracked non-sidereally 
to follow the motion of the Centaur. Field shown has North to the bottom, East to the right, 
and is 130 arcsec wide. The primary nucleus is marked with a bar. 
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Fig. 9. — The B-V color index is plotted against the V-R index, using data from Table [6] for 
the active Centaurs (blue circles) and from Peixinho et al. 2003 for the inactive Centaurs 
(red circles). Objects present in both samples are connected by straight lines. The color of 
the Sun is marked. 
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Fig. 10. — Average B — I color index vs. projected radial distance from the photocenter for 
five active Centaurs. General trend toward bluer colors at larger radii is apparent. 
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Fig. 11. — Orbital semimajor axis vs. eccentricity for active (blue) and inactive (red) Cen- 
taurs, and for Jupiter Family Comets (green). Labeled curves show the loci of orbits having 
perihelia equal to the semimajor axes of the four giant planets. Vertical dashed lines show 
the semimajor axes of Jupiter and Neptune, bounding the region of the Centaurs. 
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100 




Perihelion Distance, q [AU] 

Fig. 12. — The cumulative distributions of tlie perilielion distances of, from left to right, the 
active Centaurs in our sample, the inactive Centaurs in our sample and all known Centaurs. 
The median perihelion distances for each group are marked. 
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Fig. 13. — Derived dust mass- loss rates vs. heliocentric distance for the Centaurs in this 
sample. Circles represent measurements of the active Centaurs while inverted triangles 
represent upper limits to dust on Centaurs appearing inactive. The solid lines show power- 
laws of slope n = 2, 3 and 4, as marked, for reference (see Equation ([T])). 
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Fig. 14. — (left axis) Specific mass-loss rate [kg m ^ s ^] as a function of the heliocentric 



distance [AU] computed from Equation |8} Blue curves apply to crystalline water ice, green 
curves to CO2 ice and red curves to CO ice. Solid and dashed lines for each ice refer to 



the high and low temperature limiting described in Section 14.21 The radii of the 

orbits of the giant planets are marked for reference at the bottom of the plot. The perihelion 
distances of the active Centaurs are marked at the top with black circles. 



- 52 - 



1000 



100 — 



3^ 








5 10 15 

Heliocentric Distance [AU] 



20 



Fig. 15. — The black line marked Pk shows the variation of the Keplerian period with 
semimajor axis. Blue lines show the crystallization timescales for the lowest and highest 



likely equilibrium temperatures, as described in Section 4.2 The crystallization and orbital 
periods are equal at the points labeled A and B, having heliocentric distances of ~6.5 AU 



and 14 AU, respectively, computed from Equation (10). Left-pointing arrows indicate that 



crystallization is expected at all smaller heliocentric distances. The histogram shows the 
distribution of the perihelion distances of the active Centaurs. 



